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ABSTRACT 

Following unstable ignition of carbon, but prior to explosion, a white dwarf (WD) in a Type la supernova 
(SN la) undergoes a simmering phase. During this time, a central convective region grows and encompasses 
^ IMq of the WD over a timescale of ~ 10^ yrs, which sets the thermal and turbulent profile for the subsequent 
explosion. We study this time-dependent convection and summarize some of the key features that differ from 
the traditional, steady-state case. We show that the long conductive timescale above the convective zone and 
the extraction of energy to heat the WD core leads to a decrease of the convective luminosity and characteristic 
velocities near the convective zone's top boundary. In addition, differences in the composition between the 
convective core and the conductive exterior will significantly alter the location of this boundary. In this respect, 
we find the biggest effect due to complete ^^Ne sedimentation prior to carbon ignition. These effects add 
diversity to the possible WD models, which may alter the properties of the SN la explosion. 
Subject headings: convection — supernovae: general — white dwarfs 



1. INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmic dis- 
tance indicators has focused attention to the study of white 
dwarf (WD) explosions. Of particular importance is deter- 
mining the parameters that dictate the observed diversity of 
SNe la. Recent modeling demonstrates that the variation 
along the width-luminosity relation (Phillips et al. 1999) may 
be explained by larg e variations in the abundance of stable 
iron g roup elements dKasen & Wooslev i i2007; Wooslev et al] 
l2007h with the dominant cause for diversity likely resid- 
ing in the explosion mechanism (i Mazzali et al.ll20(>7h . An- 
other important vari able is the metallicity of the WD core 
dTimmes et al.ll2003h . 

It is critical to explore the initial conditions that may add to 
the diversity. Simulations of turbulent thermonuclear flames 
in WDs have demonstrated that the composition and energy of 
ejecta depend sensitively on the competition between flame 
propagation, instabilities driven by turbulence, and expan- 
sion of the WD ( Hillebrandt & Niemever 2000, and refer- 
ences therein). These simulations should therefore depend 
on the thermal and turbulent state of th e WD set by the pre- 
explosive convective sim mering phase dNomoto et al.l IT984t 
IWooslev & Weaverf ll986). 

During the early stages of simmering many studies have 
focused on the convective Ure a process dP aczvhski 1972 
BruennlfT97l ICouch & Arnetti IT971 Ilben[il978a.h. ,1982 



Barkat & Wheeled Il990t iMochkovitchI Il996t ISteinetal l 



1999; Bisnovatvi-KoganI 1200 U iLesaffre et al.1 l2005t 
Stein & Wheeler 2006). This occurs when nuclei re- 
peatedly electron capture and beta decay as they are carried 
by convection back and forth across the electron capture 
threshold density (the "Urea shell"). Although the energy 
loss from this process is most likely not great enough to 
cause global cooling, it ca n still have a signi ficant effect 
on the convective motions dLesaffre et al.ll2005h . Once the 
central temperature has grown above w (5 - 6) x 10** K 
(which corresponds to 10^ s before the burning wave 
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begins), there is no lon ger time for electron captures on ^^Na 
dPiro & Bildstenll2007h . and the convective Urea process will 
cease. During the last ~ 10^ s any compositional gradients 
are mixed homogeneously by subsequent convection. 

An additional place where simmering is important is 
for understanding the conditions within the WD imme- 
diatel y prior to the e xp losion (iGarcfa-Senz & Wooslevl 
1991 iHoflich & SteinI l200l IWooslev et all l2004t 
Wunsch & Wooslevl l2004l) . The properties of the tem- 



perature fluctuations present in the convection set the size 
and distribution of the ignition points, which are crucial for 
determining the s uccess of the subsequent burning wave (see 
iRopke et al.ll2006i and references therein). The interaction of 
convecti on with rotation sets the morphology of convective 
motions (Kuhlen et al. 2006) as well as the overall rotation 
profile of the WD (Piro 2008). 

The last way simmering has gained attention is in its 
abihty to enhance the neutron abundance in the WD core 
dPiro & BildstenI l2007t IChamulak et~aD l2007h . This hap- 
pens primarily via the reaction chain '^C(/?,7)'^N(e",i/e)'^C, 
where the protons are leftover from '^C burning. Depend- 
ing on the amount of carbon that is consumed before burning 
becomes dynamical, as well as the density at which it takes 
place, this neutronization enhancement could very well be 
large enough to mask any trend expected with metallicity in 
environments that have roughly sub-solar metallicity. 

In this present work we focus on the general properties of 
the simmering convection, with the aim of identifying char- 
acteristics that may introduce diversity to the SN la progeni- 
tors. We begin in ^by presenting the main features of our 
models. We illustrate how time-dependent convection in the 
simmering phase differs from the familiar case of steady-state 
convection. In this new picture, the convective flux decreases 
outside the central heating zone due to both the heating of 
new material as the convective region grows and the inability 
to transfer significant energy to the conductive exterior In ^ 
we explore the location of the top of the convective zone. We 
point out that degeneracy effects enhance the response of the 
boundary location to changes in composition. We conclude 
with a summary of our results and a discussion of future work 
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2. LUMINOSITY AND CHARACTERISTIC VELOCITIES 
FOR EXPANDING CONVECTION 

We begin by summarizing the main features of our simmer- 
ing models. (For further details, the interested reader should 
refer to Woosley et al. 2004; Lesaffre et al. 2006; Piro & Bild- 
sten 2007; Piro 2008.) Unstable ignition of '^C occurs when 
the heating from carbon fusion beats neutrino cooling. The 
central temperature then rises and a convective zone grows 
outward, eventually encompassing ~ IM© of the WD after 
~ 10^ yrs. As the central temperature, Tc, increases, carbon 
burning becomes more vigorous and the heating timescale, 
th = ((ilnTf /(if)"', gets shorter. This timescale in general de- 
pends on the size of the region responding to the rising central 
temperature at the core. 

Simmering ends and a burning wave commences once f/, < 
fconv, where fconv is eddy overturn timescale. At these late 
times, individual ed dies may experience significan t heating 
during their transit ( IGarcfa-Senz & WooslevI IT995I) . so that 
the temperature profile is no longer an adiabat and the en- 
tire convective core does not respond to the increasing central 
temperature. (In contrast, we show below that during the ma- 
jority of the time th depends on the heat capacity of the entire 
convective mass.) This makes it difficult to exactly calculate 
the precise rn oment when simmering ends. For this reason 
iLesaffre et al.l (|2006) explore fconv = cth, where a < 1 pa- 
rameterizes this uncertainty. Since we are only roughly con- 
cerned with resolving the end of the simmering phase, we 
take fconv ~ f/i ~ CpTc/e, where Cp is the specific heat capac- 
ity at constant pressure, e is the heating from carbon burn- 
ing, and all these quantities are evaluated at the WD center. 
This estimates that simmering should end when f/, w 7 s at 
a central temperature and density of T^. « 7.8 x 10^ K and 
2.6 X 10^ g cm"^, which is roughly in a greement with 



the results presented by I Wooslev et al.l ( 2004 ) using the Ke- 
pler stellar evolution code dWeaver et al.llI978h . 



We follow the simmering phase by calculating a series of 
hydrostatic WD models, each w ith a different central tem- 
peratu re, but at a fixed mass (see iPiro & Bildstenll2007l: iPirol 
|2008|) . For simplicity we ignore the convective Urea process 
since our focus is on the last ^ 10^ s. The timescale for ther- 
mal conduction across the WD is fth = Kc/R^ ^10^ yrs, where 
Kc is the conductivity and R is the radius, which is much 
longer the timescale over which heating is occurring. There- 
fore the convection efficiently mixes entropy and the convec- 
tive region nearly follows an adiabat out from the WD center. 
Outside the convective zone, we assume the WD is isothermal 
with a temperature T/. 

2.1. Convective Luminosity 

To understand how time-dependent convection is different 
than from that normally found in steady-state convection, we 
focus on the time-dependent entropy equation 

ar_ die 



(1) 



where Lc is the convective luminosity. This equation omits the 
work required to expand the WD as the heating takes place, 
which is a significant amount of energy and thus requires 
some discussion. For each convective model we compared the 
total change in WD binding energy to the total change in in- 
ternal energy of the electrons (which are primarily degenerate 



and relativistic), includin g the ion-electron Coulomb interac- 
tion energy (according to lChabrier & Potekhinlll998h . These 
two quantities are equal to the numerical accuracy of our inte- 
grations, which demonstrates that all of the work required to 
expand the WD comes from changes in the internal energy of 
the electrons. Thus, the entropy created from nuclear burning 
all goes into convective motions or the internal thermal en- 
ergy, and we are justified in omitting the binding energy and 
electron internal energy terms from equation ([T]). 

The temperature profile in the convective zone follows an 
adiabat with a power law index n = {dlnT /dlnP).^^. The 
time derivative of the temperature at a given pressure can be 
expressed as 

dT/dt = (P/PcT [dTc/dt + Tc\n{P/Pr)dn/dt] 

^(P/PcTdTjdt, (2) 

where Pc is the central pressure, and for simplicity we are 
assuming that it does not change appreciably in time. From 
this we see that the timescale for the temperature change at 
any pressure is set by the central temperature change 



dlnT _ dlnT, _ 1 
dt ~ dt ^Ji, 



(3) 



Therefore there is a well-defined, global heating timescale, 
t/,^. To account for the changing central pressure in a more rig- 
orous calculation, we must take partial derivatives at constant 
mass coordinate, M^. This can be performed by inverting the 
empirically found MdTi, Tc) relation presented in Piro (2008; 
or see eq. ifTTl below), where T/ is the nearly isothermal tem- 
perature of the non-convective, conductive region. The result 
is 



r(7;,M,) = 0.837; 
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where /i^, is the mean molecular weight per electron. This can 
be used to find (dlnT /dt)M, = dlnTc/dt, which confirms our 
conclusion that f/, is the same at any depth within the convec- 
tive zone. 

Multiplying equation ^ by dMr = Anr^pdr and integrat- 
ing, 

c„—dMr= / edM,.-Lc(Mr)+L,{0)- (5) 

/() * Jo 

We pull f/, outside of the left-hand integral to find 

^' ^dlnT 1 r^' EM) 

CpT—- — dM,-= — / CpTdMr = , (6) 

(Jt th Jo th 

where £'th(A/r) is the integrated thermal energy up to a mass 
coordinate Mr- We set Lc(0) = and define the nuclear lumi- 
nosity as Lnuc = / edMr, so that equation (|5]l becomes 



EM) 
th 



= L„„e(M,)-L,(M,) 



(7) 



Equation (|7J expresses that the nuclear luminosity must ei- 
ther go into thermal heating or convective motions, and it is 
valid at any Mr- It contains two unknowns, f/, and LdM,). We 
set the luminosity at the surface of the convective zone to be 
zero, Lc(Mc) = 0. This boundary is required since fth is long 

^ Note that this timescale is different than the local timescale, CpTc/e, used 
above for estimating the end of simmering. This is because, with the excep- 
tion of late times, the convective zone is well-coupled. 
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Fig. 1. — The convective luminosity, Lc and cliaracteristic velocities, Vc, 
as a function of mass coordinate, Mr- The upper (lower) lines in each panel 
are for a central temperature of = 8 X lO** K (6 X 10*^ K). Solid lines 
are the estimates for time-dependent convection using eq. ^9). The dashed 
lines are calculations that assume LAM, ) = Lnuc(Mr), which are plotted for 
comparison. 

in the non-convective regions, which prevents significant heat 
transfer. We can then solve for f/„ 

fA = £,h(M,)/Lnuc(M,), (8) 

which matches the definition of f/, that' Weinberg et aiTI (l2006h 
use in the context of type I X-ray bursts on neutron stars. We 
substitute ti, back into equation (|7]i to get the convective lumi- 
nosity 

L,(M,) = L„uc(M,) - E,h(M,)/th 

£,h(M,) Lnuc(M,) 



= Lnuc(M,) 



1- 



£,h(M,) L„uc(M,) 



(9) 



For steady-state convection, LdMy) = L„ac(Mr)- The ratio 
EthiMr)/Eth{Mc) is the modification due to the growing na- 
ture of the convection and the LnuciMc)/Lnuc{Mr) term is from 
long thermal time for the conductive exterior, which forces 
Lc(M,) = 0. 

In the upper panel of Figure [T] we plot as solid lines the 
convective luminosity found using equation (|9|l. These mod- 
els all have a composition of 0.5 '^c, 0.48 ""O, and 0.02 
^^Ne by mass fraction, with a mass of 1.37 Mq and an initial 
isothermal temperature T, = 10^ K. We solve for p using the 
iPaczviiskil (119831) fi t for the equation of state, an d include the 
Coulomb energy of IChabrier & PotekhinI (Il998l) . We present 
central temperatures of = 6 x 10^ and 8 x 10^ K, which cor- 
responds to tj, = 14 hrs and 170 s, respectively. The f/, associ- 
ated with the latter case is an overestimate since, as mentioned 
above, at these late times during simmering only inner por- 
tion of the core responds to the rising central temperature (in 
effect decreasing E± in eq. ||8|). The energy generation rate 



for '^C burning is taken from|Caughlan & Fowled d 19881) with 
strong screening included from Salpeter & van Horn ( 1969). 
Also plotted in Figure [T] is the convective luminosity for 
Lc(M,.) = Lnuc(Mr), i.e., for steady-state convection (dashed 
lines). Near the center, E^iMr) is small and grows less 
quickly than Lnuc(A/, ), so LdMr) is initially « Lnuc(M, ) (see 
eq. ID). At larger Mr, Lc decreases due to the effects we have 
highlighted. 

2.2. Convective Velocities 

If we take a characteristic eddy scale lc, the thermal con- 
duction timescale across an eddy is - {lc/Kf\(f yrs ~ 
(H /Rf-Kf' yrs ^ lO'* yrs, where H is the pressure scale 
height. Since this timescal e is long, the convection is efficient 
(iHansen & Kawaleil Il994l) . Using estimates from mixing- 
length theory, the characteristic convective velocity, Vconv, is 
related to i^onv via 



QglcFc 

CpT p 



1/3 



1/3 



(10) 



where Q = -{din p / dlnT)p and g = GMr/r^ is the local grav- 
itational acceleration, and we have set lc ~ H. 

In the bottom panel of Figure [T] we have plotted Vc, set- 
ting the mixing-length to the scale height, lc = H {solid lines). 
The shape of these velocity profiles are similar to those 
iLesaffre et al.l (l2005h present in the context of studying the 
convective Urea process, which is active at much earlier times 
during the simmering. The convective velocities are different 
by as much as ^ 50% from the naive estimate of Lc = Lnuc 
{dashed lines) near the top of the convective zone. 

3. THE CONVECTIVE BOUNDARY LOCATION 

Buoyantly rising eddies ascend until their density matches 
their surroundings. The boundary between the convective and 
isothermal zones is therefore set by a neutral buoyancy condi- 
tion. In practice this means that both the pressure and density 
must be continuous. If both the convective and isothermal 
regions have the same composition, the boundary is simply 
set by when the adiabatic temperature of the convective zone 
reaches t he isotherma l temperature, 7^ (i.e., where the entropy 
matches. iHoflich & S tein 2002). If the composition is differ- 
ent, a buoyantly rising eddy will be prevented from passing 
very far into the isothermal region and truncates the size of 
the convective zone. This creates an abrupt change in temper- 
ature at the boundary of size AT. In actuality this change will 
be smoothed by overshoot and mixi ng, so that the entro py re- 
mains continuous (see, for example. iKuhlen et al.ll2006l) . The 
thermal and compositional structure at this boundary will be 
complicated, especially at the last moments of the simmering, 
when relatively few convective overturns take place for any 
given location of the convective boundary. Since the convec- 
tion is always sub-sonic, we expect the overshoot to be mod- 
est and the simplification of a sharp compositional boundary 
to be a d equa te to estimate the size of the convective zone. 

iPirol (12008) presented an empirically derived relationship 
for the convective boundary, Mc, as a function of the ratio of 
the isothermal and central temperatures, Ti/Tc, 



Mc = 1.48Mfr, I — 

Pe 



1- 



1.2^ 

Tc 



(11) 



We make the substitutions M,. Mc + AMc and Ti^Ti + AT 
to solve for the change in the convective mass due to a tem- 
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TABLE 1 

Compositional Summary for Numerical Models 



Model 


Purpose 


Zone 


Xi2 


Xi3 


Xl6 


X22 


1 


Homogeneous 


Convection 


0.5 


0.0 


0.48 


0.02 






Isothermal 


0.5 


0.0 


0.48 


0.02 


2 


Neutronization 


Convection 


0.48 


0.02 


0.48 


0.02 






Isothermal 


0.5 


0.0 


0.48 


0.02 


3 


Comp. Gradient 


Convection 


0.5 


0.0 


0.5 


0.0 






Isothermal 


0.7 


0.0 


0.3 


0.0 


4 


Sedimentation 


Convection 


0.487 


0.0 


0.487 


0.026 






Isothermal 


0.5 


0.0 


0.5 


0.0 



perature discontinuity, 

AMc=-O.18M0 



AT 



(12) 



where we have scaled to a temperature ratio Tc/Ti = 8, as is 
appropriate for near the end of the simmering phase. 

The temperature discontinuity allows a small conductive 
wave to propagate out from the top boundary of the convective 
zone. Since the thermal conduction timescale (fth ^ 10* yrs) 
is long in comparison to the heating timescale (f^ ^ 10 s- 
lO-' yrs), this wave can only travel a distance ^ {th/hb)^I^H <C 
H before the growing convective zone overtakes it. For this 
reason, we ignore this detail in the following calculations. 

In the next sections we study the change in the convective 
boundary analytically. We discuss two main ways in which 
compositional discontinuities can be important. These are 
changes in the neutron abundance and changes in the compo- 
sition, whic h aff ects the Coulomb corrections to the equation 
of state. In ^3.3l we compare the models summarized in Table 
[Unumerically to confirm these analytic results. 

3.1. Neutron Abundance Discontinuity 

We first study differences in neutron abundance between 
the convective and isothermal zones. The neutronization is 
typically expressed as 



1 v^Z; 



(13) 



where A, and Z,- are the nucleon number and charge of species 
/ with mass fraction X,. The initial metallicity of the SN la 
progenitor is determined by the isotope -^Ne, which has two 
additional neutrons (Timmes et al. 2003). A mass fraction X22 
of ^^Ne decreases by an amount AF^, = 2X22/22 « 1.8 x 
10~^X22/0.02. A large enhancement of ^^Ne could be present 
in the conve ctive core if substantial gravitational separation 
has occurred ("Bildsten & Hall'2001l lDelove & Bildstenl2002l; 
lOar cfa-Berro et al. 2007). Neutron enhancement in the con- 
vective zone can also occur from electron captures during 
the simmer ing, which de creases Ye by an amount AF e ^ 
lO""*- 10"^ dPiro & Bildst en 2007; Chamulak et aLll2007l) . 

We first consider the simplest equation of state where rel- 
ativistic degenerate electrons dominate the press ure with cor- 
rections from the ideal gas of ions (in 33.21 we consider 
Coulomb corrections). 



Me 



(14) 



where K = 1.231 x 10'^ cgs and /x,- is the mean molecular 
weight per ion. 



We set the temperature to 7) + AT and mean molecular 
weight per electron to /i^, + A/i^ at the top of the convection 
zone, where T, and fig are the values of these quantities within 
the isothermal zone. We assume that /i, does not change ap- 
preciably. Note that enhanced neutronization in the convec- 
tion implies A^^ > 0. We set both P and p to be continuous 
at the convective boundary, and expand to first order. 



Ar _4/i,-m„ A^, 



-K 



4/3 /; 



(15) 



Ti 3 k^T 

Recognizing that Ef = 4-K(p/ figy^^nip, this can be written 
more conveniently as 

AT^ 1 /Xi_E^A^^ 1 /i/_E^A}; ^^^^ 

Ti 3 p.e ksTi pLe 3 p.e k^Ti Ye 

where we have defined AFp to be positive (i.e., pe + = 
Ye-AYe). Setting = 1.9 MeV(2/^^ey/^pl'\ where pg = 
p/ 10^ g cm"^, this is rewritten as 



where 7^ g = 7//10^ K and /i, = 13.7 for a plasma with X12 = 
Xi(, = 0.5. The large prefactor demonstrates that degeneracy 
greatly enhances the small changes in /ig on the temperature 
discontinuity. This is because the ions provide only a small 
contribution to the pressure, so that a large temperature jump 
is needed to offset a small change in density. Using equation 
([T2I1, we find 



1/3 



4/3 



AK 



AMe = -0.09 Mg 



Ps 



1/3 



7;,8/8 



f— ) 

V13.7/ 



10/3 



10-3 



,(18) 



where T^s = Tc/lO** K. It is interesting to note that although 
Ti affects the absolute location of Mc (see eq. ITTI ). AMc is in 
fact independent of 7]. 

3.2. Mean Molecular Weight per Ion and 
Coulomb Correction Discontinuities 

We next consider the temperature difference from changes 
in the composition of the nuclei. This is of relevance be- 
cause evolution of the progenitor star during the helium 
burning stage enhan ces '^C versus ^^O with larger radii 
("Stran iero etani2003h . 

First we present the effect of an increase of the mean molec- 
ular weight in the convection zone, A^,, using the equation of 
state given by equation ( fT4l ). Setting the density and pressure 
continuous across the convective boundary, 

Ar Au; 

Ti Pi 

As a concrete example, we consider the composition in Model 
3 from Table [T] which gives Apij pi « 0.06. This implies a 
change in convective mass of 



AMc = -1.1 X 10"^ Mq 



1/3 



0.06 ' 



(20) 



7i.8/8 \Pe 

which is negligible. 

However, compositional changes also alter the Coulomb 
contributions to the internal energy. For a multi-component 
plasma, the strength of this effect is measured via the 
Coulomb parameter 



(zV 

aksT 



p^{Z^ fU.7 
T, 48 



:21.1 



1/3 



(21) 
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Fig. 2. — Temperature profiles for the compositional models summarized in 
Table[T] For each model we consider a central temperature of Tr = 6 X 10** K 
(dashed lines) and = 8 X 10^ K (solid lines). All models have an isothermal 
temperature of 7^ = 10* K and a mass of 1 .37 Mq . The convective boundary 
is at the temperature break. 



dSchatz et al.lll999h where a is the mean ion separation de- 
fined as AiTw' p/{2>^imp) = 1 and 



(22) 



We ta ke the fitting function found by IChabrier & PotekhinI 
(Il998h in the limit T > 1 to model Coulomb effects. 



fie 



3 



(23) 



where Ai = 0.9052 is a fitting parameter. To isolate the 
compositional and temperature dependencies of F, we write 



r = Fo/r. Within the convection, F is larger by a fractional 
amount 

AF _ AFo AT 
— - 



To Ti 



(24) 



where AFq is strictly from changes in composition (i.e., 
changes in (Z^) and fii). Setting the pressure to be continuous 
across the convective boundary we find 



AT 
~T~ 



Ai^AFq 
3 Fo ■ 



(25) 



Due to the factor of F, a small fractional change AFq /Fq of 
merely ~ 10% (as present for Model 3) implies a mass change 
of 



1/3 



AM,. = -0.12 Mp 



(Z2) /13.7 



7^,8/8 48 



Ml 



1/3 



AFq/Fq 
0.1 



(26) 

Comparing equations ( fTSl ). ( l20b . and ( |26] | shows that both 
changes in and in Coulomb corrections can make non- 
negligible corrections to the convective boundary. Which ef- 
fect is largest depends on the specific progenitor model. 

3.3. Numerical Models 

We now calculate a few models to compare how the con- 
vective boundary changes due to compositional discontinu- 
ities. These models use the same microphysics as described in 
92.11 Each model is summarized in Table [T] and is motivated 
by a plausible progenitor scenario. Model 1 does not have a 
compositional discontinuity and is included for purposes of 
comparison. Model 2 shows the effects of neutronization by 
having a AYg « 8 x 10""^, compar a ble to the level of neutron- 
ization found bv iPiro & BildstenI (l2007h . For simplicity we 
assume that all the additional neutrons are in '^C (instead of 
heavier elements such as ^■'Ne). Model 3 considers a change 
in the mass fractions of '^C and '^O to show the effects of 
a compositional gradient. Finally, Model 4 shows the maxi- 
mum effect of sedimentation by assuming that all of the ^^Ne 
has sunk into the convective core. 

In Figure |2] and [3] we compare the four models for central 
temperatures of = 6 x 10*^ K and 8 x 10*^ K, and a WD mass 
of 1.37 Mq. In Figure|2]we set Tj = 10^ K, while in Figure[3] 
we set Ti = 2x 10^ K. While our models are not self-consistent 
in that we assume a fixed, non-evolving compositional discon- 
tinuity, they show that the location of the convective bound- 
ary can vary substantially depending on the progenitor model 
and the nuclear reactions that take place during the simmer- 
ing. Comparing Figures|2]and[3]shows that AM, can be large 
even when T, i s inc rea sed. This is roughly consistent with the 
analytics of § 93.1l and l3.2l that shows AMc is independent of 
Ti. 

4. CONCLUSION AND DISCUSSION 

We have highlighted some important properties of the con- 
vection present during the pre-explosive carbon simmering 
phase of SNe la. The convective velocities near the top bound- 
ary are decreased significantly because the convective lumi- 
nosity is extracted to heat and grow the convective zone and 
the long conductive timescale of the non-convective exterior 
enforces = at the top of the convective zone. The size of 
the convective zone can change depending on compositional 
gradients within the progenitor WD. 
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Fig . 3 .— The same as Fig. E] but with 7] = 2 X 10^ K. 

Perhaps the most severe effect depends on whether there is 
time available for sedimentation of ^^Ne prior to carbon igni- 
tion. There remains considerable uncertainty in the timescale 
for this process to occur For a 1.2M(^ WD, complete sed- 
imentation requires 9 Gyr (Delove & Bildsten 2002) with 
more massive WDs having an even shorter sedimentation 
time, but this all depends sensitively on the size of the dif- 
fusion coefficient used. The significant impact this may have 
on SNe la progenitors means that better theoretical calcula- 
tions of this diffusion coefficient is an essential concern. The 
work by Daligault (2006) is an important step in this direc- 
tion, but the multi-component plasma present in the WD may 
be a crucial detail. 



The presence of the convective core prior to explosive car- 
bon burning may have other implications for SNe la that de- 
serve a closer look. The convective motions should give rise 
to a massive core dynamo, which setting B^/Stt ~ pV^ implies 
fields on the order of B ^ 10'^ G. In addition, the motion of 
these convective eddies may stochastically excite waves (both 
g-modes and p-modes). These waves propagate away from 
the convective zone and transfer energy to shallower regions. 
Both of these processes may prove important in understand- 
ing the evolution toward explosive burning, and we plan to 
investigate them in more detail in a future study. 

Multidimensional simulations are needed to study the 
physics of this simmering phase not captured by our simple 
calculation. These simulations are very challenging in that the 
flow must be followed for many turnover times. As a result, 
im plici t Eulerian schemes (see for instance in two dimensions 
oflH5flich& Stein 2002; jstein & Wheelerf l2006h . anelastic 
codes ( Kuhlen et al. 2006) or a low Mach number formu- 
lation (Almgrenetal. 2006a b) are required. These simula- 
tions w ill be helpful in characterizing the flow topology (e.g., 
iKuhlen et al. 2006) and the number and spatial distribution of 
ignition points. 

The simmering phase helps set the nature of the eventua l 
run away and subsequen t explo sion. IWooslev et al.l (|2004|) 
and IWunsch & WoosleV (2004) demonstrate that the sum- 
mering is paramount for setting the initial ignition points 
for explosive burning. The large scale structure of the con- 
vective turbulence m ay also be impor tant for flame propa- 
gation. For instance, iHoflich & SteinI ([2002) argue that the 
initial velocities of the burning fronts are determined by 
the background convective motions and not by the laminar 
flame speed or Rayleigh-Taylor instabilities. This suggests 
that the initial explosive burning is likely to be off-c entered 
dGarcia-Senz & WooslevI 119951: IWooslev et al J |2004|) . The 
background turbulent state will also likely subsequently af- 
fect the motion of the bubbles by acting as a viscous drag 
(IZingale & Dursill2007h . 

Finally, the material above the convective core is devoid 
of this turbulence. The burning properties may change in 
an interesting manner as the flame passes into this relatively 
"quiet" region. Some have argued that a delayed detonation 
transition ( DDT) of the burning may be needed t o match ob- 
servations (jPlewaitiD IIQOI lUvMieti^ Although 
the concept of a DDT has been consider ed for some time 
( lKhokhlovll99Tl: IWooslev & Weaver" 1994). how and if it oc- 
curs isstflljiiKert^^ (Niemever& Wooslev 1997; Niemeve^ 
1 1 999t IWoosleyl l2007 ) . The position of the convective bound- 
ary may be an important detail. At the end of simmering it 
typically lies at a density of 3 x 10^ - 10^ g cm"^, which 
may become near the density usually invoked for a DDT 
(^ 10^ g cm"^) once the WD has expanded from the propa- 
gation of the deflagration wave. 
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